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TECHNICAL MEMORANDUM 1451 


ON THE STATISTICAL THEORY OF TURBUIENCE* 
By W. Heisenberg 



The interpretation of turbulence presented in the preceding paper 
by v. Weizsaecker is treated mathematically with the aid of the custom- 
ary method of Fourier analysis. The spectrum of the turbulent motion 
is derived to the smallest wave lengths, that is, into the laminar 
region; the mean pressure fluctuations and the correlation functions are 
calculated. Finally, an attempt is made to derive the constant which is 
characteristic for the energy dissipation in the statistical turbulent 
motion from the hydrodynamic equations . 

In the statistical theory of turbulence developed by G. I. Taylor 1 
and v. Karmdn,^ the irregular turbulent motion of a fluid is described 
by several characteristic functions between which simple mathematical 
relations exist: the "spectral" distribution of energy to waves of dif- 

ferent wave length, the correlations between the velocities at points 
along a prescribed displacement in space or time, and the like. The 
reports of G. I. Taylor contain detailed empirical and theoretical data 
on these functions. 

In his preceding report, v. Weizsaecker^ derived the most important 
of these functions, the spectral distribution of energy for the limiting 


*"Zur statistischen Theorie der Turbulenz." Zeitschrift ftir Pbysik, 
vol. 124, 1948, pp. 628-657. 

iTaylor , G. I.: Proc. Roy. Soc. A 151, 421 (1955); 156, 507 ( 1956 ); 

164, 15 (1958); l64, 476 ( 1958 ). 

2Karmj£n, Th. v. : Journ. Aero. Sci. 4, 151 (1957)* 

5c. F. v. Weizsaecker: ZS f. Phys., being published. V. Weizsaecker 1 s 

paper and the present treatise have been written in close collaboration 
during the time of cur stay in England in 1945. Only after the articles 
had been finished, Mr. G. I. Taylor kindly told us (spring 1946) that 
essential ideas in these articles had been found and published already 
by Kolmogoroff : Compt. Rend. Acad. Sc. USSR 50, 501 (l94l); 52 , l 6 (l94l); 

and Onsager (Phys. Revue 68, 286, 1945). Compare a report by G. K. 

Batchelor at the VI. Internat. Kongress F. angew. Mee hanik , Paris 1946. 
Approximately at the same time we learned, furthermore, about a paper 
by Prandtl and Wieghardt which contains similar concepts and has mean- 
while been published in the Gottingen Academy reports (Nachrichten der 
Akademie der Wissenschaften in GSttingen, Math. -Phys ikal KLasse of the 
year 1945, p. 6). The present paper may therefore be regarded only as 
a supplement and completion of these earlier investigations. 
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case of large Reynolds numbers on the basis of similitude considerations. 
The following sections of this paper first will translate v. Weizsaecker 's 
considerations into the accustomed terms of Fourier analysis and, with 
the aid of this translation, study the discontinuity of the spectrum for 
high frequencies due to molecular viscosity. Then conclusions will he 
drawn for the correlation functions and the pressure fluctuations, and 
finally, a derivation of the fundamental constant of the energy dissipa- 
tion will he attempted. 


1. Representation of v. Weizsaecker 1 s Considerations 
in the Terms of Fourier Components 

At sufficiently high Reynolds numbers the energy dissipation for 
the turbulent motion takes place in such a manner that the large turbu- 
lence elements lose energy due to the fact that, for them, the energy 
and momentum, transfer by small turbulence elements has the effect of an 
additional viscosity (cf. for instance Prandtl4). 

Under steady-state conditions, energy thus is continuously trans- 
ferred from larger to smaller turbulence elements, with the spectral 
region of a certain wave length always receiving from larger waves as 
much energy as it gives off to smaller waves. For maintenance of this 
equilibrium a certain energy distribution is necessary which, when the 
molecular friction is neglected, is represented according to v. Weizsaecker 
by the law 





F(k) -v k-5/5 

(k = — signifies the wave number, v n = J v^ or measure of 

\ wave length u V 

the mean velocity. 5 ) This spectrum F(k) ~ k"5/5 is bounded on two sides. 
For small wave numbers, that is, large wave lengths, for one reason or 
another, it will no longer be possible to regard the flow’ as isotropically 
turbulent, for the largest turbulence elements are governed by the geom- 
etry of the devices which generate the turbulence. This end of the spec- 
trum for small k therefore cannot in any case be the subject of a purely 
statistical theory. For large k, in. contrast, the spectrum is bounded 

^L. Prandtl, StrSmungs lehre (Flow theory). Vieweg, Braunschweig, 

3rd edition 194-2, p. 105 ff. 

5our Vq differs from the Vq in v. Weizsaecker 's report by a 
numerical factor of the order of magnitude of 1. 
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"by the molecular viscosity. For large k, finally, the molec ular vis- 
cosity will "become larger than the apparent turbulent viscosity, and 
the spectrum then will drop off very rapidly. 

For the calculations, we shall use the following notation: the 

velocity v in a normalized volume V is to he expanded into a Fourier 
series 


X = ^ Xk e ’^~ r (^x “ n x* • • • i &x> n y* n z integers^ (2) 
k 

Therein vjj. = v*_ k 811(1 number of the "natural vibrations" between 
— — 2 

k and k + Ak: is given^ by ^2^ — Then one obtains 

(aO 5 





F(k)dk 


thus 


F(k) = (2rt)" 2 k 2 y 



From div v = 0 there follows 


( 3 ) 


(*> 



( 5 ) 


Let us call the coefficient of the molecular- viscosity, p. The 
mean energy loss then is, because of div v = 0 with the assumption 
that the bounding surfaces are at rest: 


^This method is somewhat less apparent, but mathematically more 
convenient than the customary expansion with respect to sin and cos. 
It amounts formally to the limiting condition that v, Sv/Sx, ... at a 
bounding surface of the volume V are to have the same values as at the 
opposite bounding surface. ' 
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S = |i(rotv) 2 


(6) 




= n jTF(k)2k 2 dk 


(7) 


If the spectrum. obeys in a large region the law F(k) ~ k“5/3 } the 
total energy is determined hy the largest turbulence elements. We may 


assume for instance that the law k“5/3 is valid down to a mi nimu m wave 
number k@; for smaller k we shall assume F^) = 0* Then 


Thus 


and 



poo ~ 00 

2 / F(k)dk - 2 / <3k = 

J 0 J kQ k 5/5 


C= \ ^ 


F(k) = 


v 0 2 k Q g/3 
5 k 5/3 


for k > 


>*0 



(9) 


2 1 = vo 2 k 0 2 /3(2 J r)3 
1 6jrVk i;L / 5 


( 10 ) 


V. Weizsaecker considers the energy loss of that portion of 

the total spectrum the wave numbers of which lie below k. For these 
turbulence elements, the turbulence elements of smaller wave length 
(<2rt/k) have the effect of an additional viscosity. One may therefore 
generally write 


Sk 


(m. + J F(k') 2k' 2 dk' 


(H) 


where is to designate the additional turbulent viscosity; it i§ 

produced by the cumulative action of all turbulence elements with wave 
lengths <2n/k. With respect to dimensions, is according to 

Prandtl the product of density, mixing length, and velocity where the 
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mixing length will "be comparable to the diameter of the turbulence ele- 
ments in question whereas the velocity of the turbulence elements is 

given for instance by v 0 (ko/k)l/3. With reference to v. Weizsaecker 1 s 
report one will therefore put 


% = K P 



( 12 ) 


(k is a numerical factor) 

The expression below the integration sign is essentially deter mine d by 
dimensional considerations; but one could of course -tmaginp that, for 
instance, the waves k' in the proximity of k enter into the integral 
with somewhat different weights than the waves with large k values; 
that is, the integrand could depend, in addition, on the dimensionless 
number k’/l* However, because of the homogeneous form of the spectrum 

F(k) ~ k~5/3, one may include all these uncertainties in the numerical 
factor k and give to the integral arbitrarily the exact form ( 12 ). 

This method is unobjectionable in the region of thd k - 5/3_iaw but 
becomes inaccurate at the ends of the region where the geometry or the 
molecular friction modifies the spectrum. But even at the latter limit 
( 12 ) will still be a good approximation which at least qualitatively 
correctly represents the effect of friction. 

The constant k in equation (12) must be exactly determined by 
the hydrodynamic equations; it has the same numerical value in n.l 1 cases 
where one may speak of statistical isotropic turbulence, and does not 
depend in any way on the geometry of the flow. The theoretical deter- 
mination of this important number will be attempted in section 5. There, 
it will also be shown that the turbulent energy dissipation actually may 
be written as a double integral of the type 

pk P oo 

/ dk*.../ dk' 1 ..., with the integrand signifying the energy 

J 0 J k 

transferred from k’ into k' ' per unit time (equation ( 89 )). This 
integral is more complicated than the simplifying expression ( 13 ) which 
results from (ll) and ( 12 ); however, for the following considerations (ll) 
and (12) may be regarded as sufficient approximations. For one thus 

obtains 


St =(h + pk J dk"y F(k') 2 k' 2 dk' 


( 13 ) 


The decisive step of v. Weizsaecker 1 s consideration is the state- 
ment that this expression for k » kQ must be independent of k: 
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S]£ = S = const (for k » ko) 


( 1*0 


"because the total energy lies, almost entirely, in the long-wave region 
of the spectrum, and the energy "transport" thus must become independent 
of k. ___ 

Equation (13) may therefore he interpreted as the_determining equa- 
tion of the turbulent spectrum F(k) which must yield for the region of 

large Reynolds numbers the k“5/3-law, and for still larger k values 
the fading of the spectrum due to the molecular viscosity. 


2. Shape of the Spectrum in the Region of the 
Smallest Turbulence Elements 

We put first p/p = v and differentiate (13) with respect to k. 
Then there results 

£ + J dk' ^^^ F(k) k 2 = y&jl J F(k' )k' 2 dk' ( 15 ) 


Then we define new variables x and w by the equations 

k 


x = 


= lg jg, F(k) = F(ko)e- w ; w = w(x) 


( 16 ) 


Therewith (15) is transformed into 



(IT) 


v vk 0 \P5 , ... 

The constant h /— = — — — (cf. equation (9)) is essentially the 
K y Fo Kv 0 

reciprocal Reynolds number of the total flow and therefore always very 
small; if the Reynolds number itself were small, the flow could not be 
turbulent at all. By repeated differentiation there originates from (17) 


1 dw\[ v | p - 0 
\2~ 2 dx/U 1/F o 


w 

2 


~dx = 2e 


w x 
2 "2 


( 18 ) 
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/: 


In this equation., one may approxi m ately evaluate the integral 
w+x 

“ A 

dxe hy expanding w(x) in the proximity of x: w(x-jJ = 


v(x) + (x L -x) g + 


and breaking off after the second term. Since 


the exponential function rapidly decays, one thus obtains a good approx- 
imation 


1: 


_ v+x 
dxe ^ 


w+x 

2 


i + 4* 

dx 


( 19 ) 


By substitution into (18) there finally results 

^0 v+x 


(7 aw\/_L 

V " to A 2K v F o 


1 + 4* 


= 2 


(20) 


dx/ 


One recognizes from (20) immediately the variation of the spectrum. For 
not too large x and w the first term in the sum may be neglected and 
one obtains 


(j - S) ■ 4 + 1 } 1- £ - S» w 

as must be the case according to v. Vfeizsaecker 's theory. For large x 
and w, in contrast, the first term predominates; therefore one must 
then have 

= 7 F ~ const k -7 (22) 


In the region of the smallest turbulence elements the spectrum therefore 
decays very rapidly, namely with the seventh power of the wave number. 

Only in the transitional region from (2l) to (22) are numerical 
calculations necessary in order to determine the solution of (20). Since 
one may put for smaller x, that is, for the region 


1 


« x « 1 lg Q- 
b v 



w 



X 


(23) 
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(therewith not only (18)., but also (17) then is satisfied with sufficient 
accuracy), one may calculate, progressing from, point to point, dw/dx 
according to (20) from w, and therewith derive w for higher x. It 
is sufficient to perform the numerical calculation for a particular 

2k / Fn 

large value of the constant, for instance — ,/r-a. = a. -For another 

v VhQ 

value b one may then obtain w by a simple similitude transformation: 


wt(x) = W a (x + | lg l) - \ Xg| 


(24) 


as one recognizes by substitution into (20) and (23). 

Figure 1 shows the result of the numerical calculation for ■ 

2k / Fo 

~v~ JlSQ ~ a '' :)0 ' u ^ The numerical integration shows that for this 

value in the region 

x > 5 

w(x) « 7x - 21.85 

More generally, one obtains therefore in this region of the k“ 7 - law: 


(25) 


that is. 


w(x) = 7x + 3-0 - klgk I®) 

v V -*^0/ 

F 5, ,Vk^\ 7 

F(k) = 0.0496 (Si) 

ko 2 ' v / / 


(26) 


(27) 


A serviceable interpolation formula which is correct in the two 
limiting cases and does not result in any large errors in the transitional 
region, either, reads: - 


/ko\5/3 
F(k) = F 0 (^j 


/ k \ 8 / 5 
1 + F" 

\ k b/ 


(28) 


If one defines Lq = ft/ 3sq as the "diameter of the largest turbulence 
elements" and introduces as the Reynolds number of the total flow 


R 0 = 


PVqIo 


(29) 
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one obtains according to (9), (27), and (29) 

k s = 0.l6ko(R OK ) 3 / 4 (30) 

One may denote Lg = 3t/k s as the "diameter of the smallest turbulence 
elements" and obtains 


L s = 6.25io(Ro (51) 

By (9), (28), and (30) the form of the spectrum is determined for 
the entire k-region. For the actual flows, of course, the shape of 
the spectrum will be different for small k-values (k~ko) since there 
the geometry of the tests plays a role, for instance, the shape of the 
grids by means of which the turbulence is produced. In order to be able 
to sensibly carry out the comparison with experiments, one will then 
introduce a quantity ko in. such a manner that, for instance in the 

domain of the k - ^^-law (thus for k s » k » ko), the formula 


F(k) 


3koW 


(52) 


becomes correct. The quantity kQ thus defined then does not give any 

direct statement regarding the variation of the spectrum for the smallest 
k-values. Generally, though, the spectrum will greatly deviate from the 
-5/3 

k -law in the region k ~kQ. 


For k » k B , too, the spectrum will not unlimitedly retain the form 

k“Z. The well-known investigations of Burgers/ make it very probable that 
at sufficiently small Reynolds numbers finally no turbulent motions what- 
ever exist. On the other hand, the k“/-law shows such rapid decay that 
the region k » k s is practically insignif leant . A somewhat larger 
error will arise, particularly in the transitional region, due to the 
inaccuracy of equation (13) itself; but it is probably not worth while 
to apply already at this point the much more complicated equation of 
section 5 to the problem stated here. The correct equations would at 
any rate lead to somewhat different numerical factors in (27), (30), and 
(31). 

/Burgers, J. M. ; Verh. d. Kgl. Nied. Akad. d Wiss. 17, Nr. 2, 1 
( 1939 ) > 18, Nr. 1, 1 (1940). 
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For the comparison with experiment one needs the energy distribution 
with wave numb er for a certain direction, for instance with since 

the spectra have been measured experimentally by Simmons 8 and Dryden9 
by means of the fluctuations of the velocity with time in an airstream 
which is guided past the measuring point at a constant velocity U which 
is large relative to v Q . Also, this spectrum has a different form 
according to whether the Fourier expansion of y x or Vy is concerned. 

Exper ime ntally, first the spectrum for v x is required; however, we 
shall also derive the spectrum for v y since it will be necessary 

later on, in the calculation of the Correlation functions. Since 
according to equation (5) (v^k) = Oj one obtains 



The spectrum of v^. in k^. which we shall designate by F x (k x ) 
becomes therefore 

‘if $ (i 2 - 

In a similar manner, the spectrum for Vy. becomes 


(34) 


= t/J§( k2 + k * 2 ) F(k) 

From (34) and (35) there follows for 
( a ) kQ « k « kg : 


(35) 


F(k) - F 0 (^) 5/5 thus F x (k x) - ^fe) 5/3 and Fyfk*) - 


(36) 


^Simmons, U. Salter: Proc. Roy. Soc. A 165, 73, (193 8). 

^Dryden, Schubauer, Mock u. Skramstad: Nation. Adv. Comm. Aero., 

Nr. 581, 1938; Dryden, H. L. : Proc. V. Intern. Congr. f. Applied. Mech. 

Cambridge, Mass., p. 3^2, 1938. 
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(b) k » k s : 

*.5/3*16/3 

F(k) - Fq ^ s 

*i(M - Ss *o 5 / v- 6 / v 7 

Fy(k x ) = Jfi (37) 


As a serviceable interpolation formula (which, however, is somewhat 
less accurate than ( 28 ) in the transitional region), one may again put 


with 


and 


with 


F /v\ - E 0 /V \ 5/3 

*x(rxj llolk x J 


1 + 


k x \8/3 


iKr / 


-2 


k^ = 0 . 6 * 1-5 k s 


- S^f 5 


1 + 


^ w 5 




-2 


= 0.793 


(38) 

(39) 


(to) 

(*U) 


Before the comparison with experiment is carried out in detail, we want 
to raise the question at what critical Reynolds numbers the transition 

from the k‘ -5/3 

-law to the k“7-law takes place, that is, - if one wants 
to express it in this manner - the transition from the really turbulent 
motion proper to the laminar motion. One may regard as the critical 
Reynolds number for this, perhaps, the expression 




wherein, according to v. Weizsaecker, v g 


and ( 31 ) one then obtains 


*s = 


10.2 

K 



(42) 

From (29), (30), 


(*3) 
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The numerical value of k will be discussed later on. At any rate, the 
transition therefore takes place at a certain numerical value of the 
Reynolds number, as was to be expected from general similitude considera- 
tions. 


In figure 2, the measurements of the spectrum F^k^.) by Simmons IQ 

are compared with the theory. The measurements in question are intensity 
measurements on an airstream which flows past the measuring point at the 
velocities U = 456 cm/sec (o), 6o8 cm/sec (X)> 1060 cm/sec (□), and 
which has been made turbulent by a grid of J.6 cm mesh width; the meas- 
urement was made 2.1 m behind the grid. The measured points of Simmons 
are plotted individually only in the right part of the figure, in the 
left half the approximate scatter of the measured points is indicated 

by a vertical line. The abscissa is k in cm - -*- (in logarithmical 
scale), the ordinate F x (k x ) , likewise logarithmically, in arbitrary 

units. If one assumes U/vq to have the same value in all three meas- 
uring series, which is confirmed by other measurements by Taylor, one 
obtains, in the case of a suitable selection of this ratio, the three 
curves plotted in the figure. Qualitatively, the Simmons data are well 
represented by the curves, particularly also the divergence of the three 
test series in the short-wave part of the spectrum. In details, however, 
there exist considerable discrepancies; one recognizes from the figure 

that the range of validity of the k - 5/5 

-law Is here so small that a 
reliable check Is not possible. The reason for this is the smallness 

of the Reynolds number Rq. For k = 1 cm"- 1 - the diameter of the turbu- 
lence elements is 3 cm, thus about half as large as the mesh width of 
the grid; in this region, the turbulence Is not yet fully Isotropic, 

therefore the k~5/5 -law cannot yet be valid. However, already for 

k = 4 cm”! the Influence of the molecular viscosity becomes noticeable, 
and the intensity drops off markedly. The related measurements by Dryden 
quoted before which extend over a large spectral region have been made 

at Reynolds numbers so small that the validity of the k - 5/5 -law can 
hardly be checked. Therefore it would be desirable that similar meas- 
urements be carried out at very much larger Reynolds numbers. For the 
ratio u/vq, one obtains from the adjustment of the theoretical curves 

to the measuring points of U/vq = 55 k if one identifies I q with the 

mesh width of the grid. This value agrees well with measurements!! of 
this ratio in similar tests if one assumes k. to be about 0.5. 


10 Simmons, U. Salter: Proc. Roy. Soc. A 165, 73 (1938). 

!lCf. G. I. Taylor, Proc. Royal Soc. A, 164, 486, (1938). 




NACA TM. 1431 


13 


Another and probably more accurate determination of k is obtained 
from the damping with time of the turbulence, already theoretically treated 
completely by Taylor (cited before). For the total energy loss per cm? 
and second: S, there results from (9)> (13)> and (l4) 


s = P* “P ▼<&> 


(44) 


For the damping with time of Vq there must thus apply 


_d_ 

dt 



= ‘ v 0 5kQ 


with the solution^ 


(^ 5 ) 


v 0 (t) 


▼o(°) 


1 + -5^- K kov 0 (0)t 


(46) 


Taylor who essentially derived this equation reports on measurements 
by Simmons in which u/v 0 (t) was ascertained as a function of t = x/u 

(x equals distance of the measuring point from the grid). From (46) 
follows 


U 

v 0 (t) 


U , \/~3 x _ U 

v 0 (0) 8 Lq v 0 (0) 


+ 


0.68k 


X 

*0 


(^7) 


If one puts Vq = u ' \p5 (u' = y v x 2 according to Taylor) and identifies 

Iq with the mesh width, there follows from Taylor's measurements k = 0.85, 
from the corresponding measurements by Dryden a somewhat smaller value. 
However, because of the uncertainty regarding the value to be inserted 
for Iq* this determination is probably still uncertain by about 50 per- 


cent. 


3- The Correlation Functions 

Taylor and von Karma, n (cited before) studied the correlations which 
exist between the velocities at two points at a given distance. The two 

12 Footnote at the time of proof correction: For this solution, 

kQ = const, is presupposed which certainly is not the case for larger 
times. The problem of damping is investigated more closely in a paper 
of the author about to be published (Proc. Roy. Soc. A.). 
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correlation functions Rj_(x) and R 2 (x) which therein play the main 
role, are defined as 


B l( x) , Zgspg) % (x) = M 


m 


with the point P 2 displaced with respect to the point P-j_ by the 
distance x in the x-direction. 

These functions are, according to Taylor, in a simple relationship 
with the spectra: 


Rl(x) 


_ 0 


dX^kxJcos k^ 


d^xFx(-^x) 


R2(x) = 


poo 

/ dkxFy^hxjcos k^x 

J 0 


(**9) 


p co 

J o ^FyC^x) 

With the aid of equations ( 3 ^) and (35), (^9) is transformed into 

r oo 

dkF(k)(sin kx - kx cos kx)k"5x"5 


R l (x) = 


p 00 

/ dkF(k) 

J n 


R 2 (x) = 


— / dkF(k)(k2x2sin k x + kx cos kx - sin kx)k“5x"5 

0 _ 


p 00 

/ dkP(k) 
J n 


(50) 


From these expressions one recognizes immediately the correctness of 
von Kerman's relationship 


R 2 = 


Ri 


+ x£5l 

2 dx 


(51) 
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The formulas (50) may "be approximately evaluated in the two limiting 

cases x « l/k s and l/k s « x « l/kQ. When x « it is advisable 

K s 

to expand the integrands with respect to powers of x. The first terms 
of the expansion then lead to the quantity 


k2 = 



(52) 


which may easily be calculated from (13), (44), and (29): 


k 2 = KEQkQ 2 


Thus one obtains: 


R x (x) = 1 - 


x 2 k 2 

10 




For x « l/k s 


/ \ x 2v2 

R 2 (x) = 1 - £JL- + . 


Taylor defined a length A by the equation 


(55) 


(54) 



(55) 


and designated it as a measure for the magnitude of the smallest turbu- 
lence elements. According to (55)* A becomes 


A = 2.71 = 0.434l s (r oK )4 (56) 

It mus,t be emphasized that A is not identical with the quantity 
L s (equation (5l)) which we have denoted as "diameter of the smallest 

turbulence elements" and that A also depends on Iq and Vq in a 
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manner different from that of Lg. A comparison of (31) and. ( 56 ) shows 
that for sufficiently large Reynolds numbers Lg becomes «A. 

In the opposite limiting case l/k s « x « 1 /]£q one obtains 
from (50) 


Rx(x) = 1 - 0.643(kQx} 2 /5 + . . . 

(57) 

R2(x) = 1 - 0.858(kQX)^/5 + . . . 

Here these first terms of the expansion do not depend on the special form 
of the spectrum in the proximity of kQ ; only for x ~ l/kQ does the form 

of the spectrum in the proximity of ko becomes important j there, however, 

the problem may no longer be treated with purely statistical methods. 

The formulas ( 5 ^) and (57) thus give essentially a complete description 
of the correlations in so far as they may be regarded as a consequence 
of statistical isotropic turbulence. 7 

The formulas (5^0 and (57) also show clearly that the correlation 
function does not have the same form in all flows but that rather, in the 
case of variations of the parameters, the inner and outer parts of the 
function undergo different similitude transformations. ThiB point has 
been stressed particularly by Taylor 15 in contrast to a different con- 
jecture of von Karmen (cited before). 

For comparison with experiment, the measurements of Rj_(x) and Rg(x) 
made by Simmons have been plotted (circles and dots, respectively) in 
figure 3; furthermore, the theoretical curves calculated according to 
the exact formula (50) are shown. Here again Lq has been identified 

with the mesh width of the grid and A has been calculated from the 
spectrum for U = 1060 cm/sec. The experimental values_agree, in fact, 
with the theoretical ones very exactly at the smaller values of x, 
actually more exactly than could have been expected in view of the uncer- 
tainty of Iq. Beginning from xk^ ~ 1 the deviation of the experimental 
points from the theoretical curves, becomes noticeable which was to be 
expected from the derivation. The variation for larger ”x- values depends 
on the behavior of the spectrum in the proximity of ko which cannot, 
in principle, be represented by our formulas/' But even for larger x- values 
the deviations from the theoretical curves remain small. 


-'-5g. I. Taylor, Jour. Aero. Sci, k, No. 8, 3H, 1957* 
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4 . The Pressure Fluctuations 

While studying the dif fusion in a turbulent airstream, Taylor (cited, 
before) has derived a relationship between the correlation function for 
the diffusion and the root mean square value of the pressure gradient. 

We shall investigate therefore also the root mean square values of the 
pressure fluctuations from the viewpoint of the theory here described. 

With reference to equation (2) one may expand the pressure into a 
Fourier series 

P= = p k* (58) 

k 

and the fundamental hydrodynamic equation 


▼ - - (zV)v - VP + ~ Av 


is transformed into 



% 


P 

k2 


2Z (Yk’^)(lk-k’^) 

k* 


( 59 ) 


(60) 


For the root mean square values of pressure and pressure gradient there 
results 



grad^p = 



(61) 


We are interested, first, in this latter mean value: 



18 


NACA TM 1431 


grad 2 p = XIXI XI 7 p (Yk'^KYk-k'kXZk 1 «£)(^-k' *k) (62) 

k k' t^ 77 ^ - - 


The superscript bar indicating approximation to the mean signifies 
here simply approximation to the mean in time. If one wants to calculate 
mean values of the type ( 62 ), one must take as a basis_some kind of 
"assumption of disorder" regarding the turbulent motion. One may start 
from the fact that the amplitudes vjj. in the course of t ime fluct uat e 
by a value given by (10) or (28), respectively, so that the time average 
v^.v_ k is given simply by the spectrum ( 28 ). The phase's of the v. , 

however, will in the course of time pass through all possible values; all 
values of the phase will occur, on the average, with the same frequency. 
If one could regard the phases pertaining to different wave numbers as 
completely independent statistically, there would, in taking the mean of 


such products of four factors 


-kl-2^3-kV 


be left only the terms in 


which every two wave numbers are equal and opposite; thus, terms of the 


type 




and these mean values could be replaced by the pro- 


ducts of the mean values of the squares of the amplitudes: 




( 63 ) 


Actually, however, statistical correlations will exist between the phases 
pertaining to different wave numbers since the waves mutually influence 
one another. In section 5 we shall attempt to estimate, such correlations 
in a simple case. In spite of the presence of the correlations, however, 
we are inclined to assume that in a sum of the kind (62) the terms of the 
type (63) make the largest contribution since their mean value is differ- 
ent from zero even in the first approximation, without any assumption 
regarding the behavior of the waves with respect to time whereas the 
other mean values attain a value different from zero only because of the 
finer fluctuations of various waves. We believe therefore that one obtains 
a serviceable approximation if one takes only the terms of the type (63) 
in (62) into consideration. Then there results 

2p2 (\ X ) ( z -k' x ") (--£ A ') 

k'k' * (k' - k") 2 


grad^p 


( 64 ) 
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In taking the mean with respect to the directions of* the vjc ” one assumes 
again all directions perpendicular to k for to he equally probable - 

it is expedient to use the relation 

(vk§.)(y_tb) = ab - (65) 

Thus one obtains 



If one sets (k'k' 1 ) = k'k' 1 one can perform the integration with respect 
to £ and finds 


= ~ y t | v k> 2 | | v k* ' 2 | k ' k ' 




1 + f S 2 + n s 4 . b 6 + * T ii| T 


For 0 < s « 1 there applies approximately 


*(s) » 8s (1 - 1 e 2 + si) 

v ' 15 \ 7 21/ 


If one transforms the sums into integrals and substitutes the spec- 
trum (28) into (67); there follows finally 



( 70 ) 
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One recognizes from (70) that the integrals converge at small k- values , 
and that one may therefore perform the integr ation from k = 0 without a 

considerable error. This shows that grad^p is actually determined by the 
behavior of the spectrum at large k, that is, by the smallest turbulence 
elements . Had we calculated p2, we would have found, on the contrary, that 
the integral diverges at small values of k, thus that its value is deter- 
mined entirely by the largest turbulence elements. Therefore, the value of 
p2 c ann ot at all be calculated according to the method. used here; for, 
first, in the case of small k, the spectrum has a form dependent on the . 
geometry; and second, it would surely by quite unjustified to consider, for 
the largest turbulence elements, only the mean values of the type (65) since 
the geometry certainly impresses definite phase relations upon the system of 
the largest eddies . 

Equation (70) now becomes 


grad 2 p = 

9 


b V /5l % 2/5 f o ^5 (j. + t 8/3J^ + ,8/5] 


(71) 


The double integral at the right was estimated, according to a gra- 
phical method, to be O.763; thus there follows finally (cf. (30)) 


grad 2 p = 0 . 17p 2 v Q i Hc 0 ij -/ 5k s 2 /5 


(72) 


= 0 . 05 p 2 v Q i| lc 0 2 ^K 


(75) 


Taylor (cited before) had expressed th e conJecture_t hat grad 2 p 
should have the same order of magnitude as p 2 v Q 2 (3v/Sx) 2 -thus the order 
of magnitude p 2 VQ^kQ 2 /^k B ^/5 # one recognizes now from (72) that grad2p 
must be considerably smaller, the more so, the larger the ratio kg/k^. 
The length Atj defined by Taylor: 



(7*0 
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that is, 


grad 2 p 




(75) 


must, therefore, for large Reynolds numbers, "become considerably larger 
than the length X of equation (55)* From (75 ) > (75) , and ( 56 ) follows 


X^ = X0.42^ 


(76) 


It is true that this result does not agree with the experimental 
findings. Taylor indicates, for a test which Simmons had performed 
adjoining similar experiments by Schubauer^, X^/X = 0.5 where one must 

assume approximately ^/rqk « 3-9 according to the test conditions. 

Thus one must raise the question whether the result ( 72 ) has perhaps been 
falsified by the fact that only the terms of the type ( 63 ) were taken into 
consideration in taking the mean. However, one can easily see that this 
may perhaps affect the numerical factor in ( 72 ) but that the dependence of 

kQ and k s , that is, the dependence of Xrj on (RqO is in no way con- 

nected with this approximation. For already equation ( 62 ) shows that on 
the right side the normalization factor Vq^Cq^/^ must appear, because 

of equation (10). After the mean has been taken, this factor is supple- 
mented by a factor of the dimension k2/5 which obviously can be at most 
of the order k Q 2/3 j it must be, because the pertaining integral with 
respect to k would diverge like k 2 /3 If the decay of the spectrum 

would not set in for k ~ k B with k~7. There would remain the possi- 

bility that only the numerical factor In ( 72 ) has been estimated as too 
low due to the consideration solely of the terms of the type ( 63 ). But 
it Is hard to imagine that the correct expression would increase by more 
than a tenfold - which would be necessary for Interpretation of the 
experiments . 


Perhaps the contradic tion ma y be cleared tip in the following manner: 
The main contribution to grad^p stems from wave numbers of the order k s , 

thus, from turbulence elements whose diameter measures a few millimeters. 
In Simmons' test, the airstream is heated by a hot wire of 20 cm length 
stretched across the wind tunnel, and then the distribution of the heated 
air is measured at a certain distance behind the hot wire. Precisely the 
smaller distances (5 to 15 cm) are decisive for the determination of X^- 

^Schubauer : Rep. Nat. Adv. Comm. Aero. Nachr. Nr. 524, 1935- 
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It suggests itself to assume that the hot wire itself produces in 
the airstream a small vortex street and additional turbulence, with the 
turbulence elements probably having a length of a few millimeters - that 
is, the wire increases precisely the intensity of the tur bulence in the 

spectral region which exerts the strongest influence on grad^p. The 
additional turbulence in the immediate proximity of the wire is probably 
much more intense than the original turbulence of the same wave range. 
However, this additional turbulence is rapidly damped, of course, and it 
is surely difficult to estimate whether this turbulence alone can explain 
the discrepancy between ( 76 ) and the empirical X-^-value. 


5 . Energy Dissipation for Normal Isotropic Turbulence 

The investigations of the preceding section are already closely con- 
nected with a basic problem of the statistical turbulence theory: namely, 

with the determination of the energy dissipation in the case of normal 
energy distribution, that is, the determination of the constant k in 
equation (12). In this problem the molecular friction may be neglected 
entirely. The fundamental hydrodynamical equations could therefore be 
presupposed in the form 


v = 



Vv = 0 


Furthermore, a suitably slngled-out partial volume of the fluid is to be 
selected as normalization volume which under certain circumstances Is 
moved simultaneously with the fluid, corresponding, for instance, to the 
mean value of the velocity with respect to the volume. We assume therefore 
that the volume moves with the velocity u. Then ( 60 ) is transformed into 




-iT~ ( 2 k*) 

k* 


— k-k' ~ k 2'' JL k-k 


, ,k) 




(77) 


For the calculation of the energy dissipation, one has to ascertain 
how the intensity |^k^| a certain natural vibration ^or perhaps 

better: the sum of such squares of amplitudes with respect to a small 

k+Ak \ 




spectral region Ak: 


varies in the course of time. One 
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recognizes from (77) that one needs for this purpose time averages with 
respect to products of the type 




(78) 


wherein k* = -k-^ - kg. Because of the statistically uniform distribution 

of the phases, these mean values would disappear if there would not exist 
statistical correlations between the phases pertaining to different k 
which stem from the mutual influence of the various waves, as has been 
explained already in section 4. In order to ascertain these correlations, 
one must somehow express in the equations the influence exerted upon a 
wave with given k by waves with a different k; one may do this for 
instance by representing one of the three amplitudes in ( 78 ) as a time 
integral over v^. and expressing v^. in turn by a sum over two other v^ 
according to (77)- Then one obtains products of four amplitudes v^. each 

of which, however, must partly be taken at different times. For such pro- 
ducts the considerations of section it- are valid according to which one 
obtains a first approximation by taking only products of the type ( 63 ) 
into consideration. Of course, one could continue the procedure in prin- 
ciple and attempt to calculate the other mean values of quadruple products 
by tracing them back to six-fold ones, etc. However, such calculations 
would probably become much too complicated^ the higher terms probably also 
would make a lesser contribution, and we shall thus be content with the 
first step. 

For these 'calculations, one will obviously need mean values of the 

type 


%(t)v_^(t + t) 


and we define therefore 


R k (t,T) 


k+£k 


^ %(* + |)s-k(t - §) 


k+Ak 


E vZ 

k ~ “ 


(79) 
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The summation over a small spectral region Ak has been included 
into the definition of R^(t^r) so that the magnitude of the normaliza- 
tion volume does not directly enter into R^ (t,T ) and that the mean is 
taken equally over all directions of k. Evidently the spectral region Ak 
must be selected wide enough that many natural vibrations of the normaliza- 
tion volume still may be accommodated in it (that is, k^AkV » 1), yet 
very small compared to k itself. These requirements are for the turbu- 
lence elements of the order of magnitude V themselves no longer compatible 
but for those, the- statistical methods cannot be applied anyway. The whole 
procedure thus can be carried out only if it is found that the large turbu- 
lence elements practically do not any more contribute to the mean values to 
be investigated. 

In order to obtain from the hydrodynamical equations information 
about the quantities Rj c (t,T), it suggests itself to examine the following 

expression: 



. ^ 5 - 

i(uk)vk 2 



(80) 




In this expression one can replace y^i by a time integral with respect 


. H 

to vv i : 



dt' + Vfc.(-T) 



+ %>(-T) 


(81) 


If T is selected sufficiently large, the correlation between 
y_k i ( t ) and Vjj.,(T) will disappear; it is therefore expedient to per- 
form, after substitution into (80), the limiting process T — >«>. If, in 
addition, one takes the mean with respect to the directions - thus elimi- 
nating the term with u - one obtains from (80) and (81): 
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(k k') 


t+J-T' 
-k" 

k' 2 



( 82 ) 


If one furthermore, as in section 4, takes into consideration the terms 
of the type ( 63 ) only, also replaces k - k* by k’ , and integrates with 
respect to the cosine of the angle between k and this vector, there 
follows : 


k+Ak 


- 0 ■ r 1 ^ - t 0 


H k r(t + T a T ‘, T'jlt~\'~^(lt g - k’ 2 ) jaac'^ + k’ 4 - | k^k' 2 ) - 

(k 2 + k^Kk 2 - k' 2 )lg|- ~ + | -| c®5) 


This equation presents the possibility of expressing the differential 

. dEi,(t,T) dR k (t,T) 

quotients with respect to time — — and — by the R^ them- 

dt dT 

selves^ when the Rjj. are known, one can, moreover, calculate the energy 

dissipation from ( 83 ), putting t = 0. For this purpose we shall assume 
that the entire turbulence phenomenon is either stationary or is damped 
very slowly so that the times during which the intensity | | noticeably 
varies axe very long compared with the fluctuation period of v^. The 


notation 


lYk* - 1 represents, therefore, the mean value over a time which is 


certainly much longer than the fluctuation period but is very much shorter 
than the damping time. 

dR v (t,T) 

The equation for — gives a measure for the fluctuations of 

dt 

the quantity R^.(t , t) as a function of t about its time mean: 


R k( T ) = RkC^^) 


(84) 
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One may surmise that these fluctuations are small in the region of small t 
which is determined by the small turbulence elements , and that they increase 
with growing t; this question will be further investigated later on. 

Before carrying out the further calculations one has to determine... how 
the partial volume V and its velocity u are to be chosen. One could 
first try to put u = 0 and to identify V with the total volume. How- 
ever, one would obtain an erroneous picture of the actual conditions: The 

decrease with time of the correlation function Rk(t,-r) = R^(t) as a 
function of r is determined in this coordinate system by the largest 
turbulence elements, and is therefore very rapid. One_ can show that the 
correlation function in this coordinate system is given with sufficient 
approximation by 


R k (r) = 1 


k 2 v2 T 2 
kv 0 T 12 

W 


n kv o T 

\/l2 


JO 


dxe x 


The calculations which have led to this expression need not be discussed 
in more detail since the expression is not further used later. The 
physical interpretation of the expression is given by the following con- 
sideration: The function R k (-r) in it decreases after a time of the 

order ; that is the time during which, for instance, precisely an 
KVq 

eddy of the wave length due to the high velocity in the largest tur- 

bulence elements, passes by the point of observation. The fact that the 
correlation function decreases after that time signifies therefore simply 
that the velocity in the largest turbulence elements is. of the order of 
magnitude v 0 , but, statistically, fluctuates about values of this order. 
This phenomenon is not connected with the disintegration of the eddy of 

the wave length £&. If the F(k) ~ k“5/3_i aw ls valid, it is rather to 

be expected, according to the similitude considerations of v. Weizsaecker, 
that the disintegration of the eddy takes place only after time intervals 

of the order 2 «Vq ’'"k ^^kg . On the other hand, however, the energy 

dissipation is connected with the disintegration of the eddies, not with 
the motion on a large scale . If one wants to describe the disintegration 
of the eddies in equations , one must move the coordinate system at the 
same time. One then must make the linear dimensions of the partial vol- 
ume V somewhat, but not very much larger than ^ and move it simul- 

k 

taneously in proportion to the mean velocity within it. We shall assume 
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experimentally that one can give for every k a volume V corresponding 

to it in such, a manner that Vk^ becomes independent of k and that for 
the thus selected, simultaneously moved volume elements, the correlation 

function R(t,T) is a universal function of the variables VQk ^'^kg' *t 
and v 0 k 2 / 5 k^/^T as is to be expected according to v. Weizsaecker 's 
similitude consideration. We shall show that the relations for Rk(t,T) 
arising from (83) actually can be satisfied by this assumption if the F k 
distributed according to the k“ 5 / 3 _iaw. 


are 


If one puts R k (t,-r) = g(£,n) wherein 


l - > 2/5 ^ /3 t n - > 2/3k o /5T V - y 


( 85 ) 


there follows from (83) 


5s( fer 1 - &X" X" a, '[4 - - , ) b (( ; + yy-y 13 ^ 13 ) - 




aji(y) JT an 1 s(l - £.1’ - l)s^(t + n f 11 ') y 2 ^.l'y 2 ^ 3 ] - 


> (86) 


.(t - fc* ♦ n^t - 


where 


f(y> ' 


2 y^l - | y 2 + y^ - 


(1 - y2)(l - y 4 )lg H 7 J 


( 87 ) 



28 


NA.CA TM 1431 


These equations actually do no longer contain the constants Vq. 

The reason is that the integral with respect to y converges for small as 
well as for large y; f(y) disappears sufficiently for y = 0 as well as 
for y — >». Therefore one may take the integral over k 1 , instead of from 
k^, simply from 0, without considerable error; moreover, the convergence 

of the integral for large values of k shows that the molecular friction 
actually is of no importance in this problem; the behavior of the spectrum 
in the region of the smallest turbulence elements is unimportant for the 
correlation functions R(t,r) and the energy dissipation for medium 
k-values . 

Before attempting the numerical solution of (86), we shall use the 
equation ( 83 ) for calculating the energy dissipation in the approximation 
here aspired to. For this purpose we put r = 0 and integrate the equa- 
tion ( 83 ) with respect to k between two arbitrary limits Kp and 



The integrand on the right side is an anti symmetrical function in k 
and k ' . If one calculates the variation with time of the total energy, 
that is, if one puts Kp = 0, K 2 = =°, there results therefore zero, as 

far as the integral on the right converges at all. That iB, the total 
energy is constant in time; this is a necessary requirement since the 
molecular friction has not been taken into consideration. However, if 
one considers the variation with time of the energy which is contained in 
the part of the spectrum lying between Kp and Kq, the integral on the 

right side may be transformed in the following manner (we shall call the 
antisymmetries! integrand simply J): 






an'(-j) 


d_ 

dt 


(89) 
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In the first of the two integrals on the right J is always posi- 
tive; in the second (-J) is always positive. From this notation there 
follows that the first integral may be interpreted as the energy which, 
per unit time, flows from smaller wave numbers (k' < Kp) into the region 

between Kp and Kg; the second integral as the energy which flows toward 
larger wave numbers (k 1 > K 2 ). If one puts, in particular, Kp = 0 and 
Kg » kQ, the second integral represents the entire energy dissipation; 

for the normal - spectrum [F(k) ~ k“ 5 / 5 ] it must prove to be independent of 
Kg. Thus one obtains from ( 83 ), (86), and ( 89 ) for the energy dissipa- 
tion the expression 


Ko 

n_2 

k" 


1^°*° Jks & Jo ayf<y) Jo lnB ( C " H 8 ((? ■ t> 2/5 ''i* 2/5 ) 


= ^ v o^o J 1 d y(- 1 ey) f (y) J' d ne 


(90) 


This expression is actually independent of Kg as it must be. Since 
the entire energy dissipation according to equation (44) is also given by 


P fC 


CTO J 


there follows 


k = 


5/5 i 


dy(-lgy)f(y) dT is^ - |^n)g ^ ( 91 ) 


From this equation k can be calculated numerically if the function 
g(£,T}) is known. 

We now turn to the treatment of the equation system (86). This sys- 
tem represents a considerable simplification compared to the initial 
equation (77) in so far as it does not any more contain any dimensional 

-3/3 

quantities, and is already derived from the equilibrium spectrum k ' . 

On the other hand, (86) also still contains statements regarding the 
fluctuations of the g(£,T]) as a function of the £ and is, for that 
reason, doubtlessly too complicated to permit rigorous solutions . One 
could attempt to completely neglect the fluctuations in a first approxi- 
mation, and to calculate with the mean values. Unfortunately, however, 
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it turns out that the contribution of the fluctuations in certain regions 
is large. One recognizes this from the second equation (86). 

For, if one puts 


g(^Tj) = g(t]) + Ag(£,Ti) (92) 


where g(iy) signifies the mean value over £ of g({;,Tj): 


s(tj) = g(^4) 


(93) 


there follows from the time mean of the second equation (86): 


8 


Poo 

d yf(y) / ^'^(t)' - Tl) + g( TJ 1 + Tl)J g^Tl'y 2 /^^ 


-If 


dyf(y) 'dtj 1 

UO 


*s(fi - - n)*((£ + 


- \> t \' + - 11 * 11 



(94) 


This relationship shows that the fluctuations Ag(£,q) cannot 
always be small. It is true that the left side of (94) vanishes for 
T) =■ 0 $ this follows from the relation 


f(^) = -y 4 / 5 f(y) 


(95) 
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which, will have to be discussed later, and signifies that the spectrum 

k"'^ is actually in equilibrium; however, for larger |q| the left side 
assumes appreciable values. Therefore, it is doubtful whether one will 
obtain a sufficient approximation if, in the transition from (82) to (83), 
one takes only the mean values of the type (63) into consideration. How- 
ever, I did not succeed in improving here the approximation or in 
obtaining more than a very crude estimation of (86) . 

One may perhaps assume for such an estimation that, for large values 
of q, the first term of the summation on the right side of (94) is much 
larger than the second. For in the first one, the integral is taken over 
the region q' ~ q which probably contributes a great deal whereas in the 
second, for large q, the factors Ag have already strongly decayed in 
the entire integration range. One may therefore attempt the assumption, 
at least for large ri, that the second term of the summation on the right 
side of (9!) may be neglected. In this approximation the time average of 
the first equation (86) then becomes, with use of ( 94 ): 


— ^ = -!■ f (fy <ay) 

dq 8 


r 


aVsU + *n 




(for q » 0) (96) 


One may utilize this equation, for instance, in such a manner that 
one assumes a plausible form for g(q) , leaving the scale in the 
q -direction undetermined at first, and determining it subsequently so 
that the equation (96) is valid as exactly as possible for large q. In 
this way one will describe the steepness of the decrease of g(q) for 
large q with some correctness, and precisely this steepness is decisive 
for the value of k . . 

In the practical execution of the calculation it is expedient to 
introduce, in the place of y and f(y), new variables 


s 



cp(s)ds = f(y)dy 


( 97 ) 


Then there applies as one recognizes from (87) (compare also ( 95 )): 


^(1) = - sc P(s) 


( 98 ) 
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This equation is based on the fact that the energy dissipation from 

the wave number ik to the wave number k coincides with the one from k 
a 

to ak, except for a factor qualified by the similitude transformation. 
Furthermore, there then applies in good approximation _ 


qp(s) = fs(l - s^) ^1 - j=s^ (for 0 < s < l) (99) 


and for larger s one can reduce cp(s) by means of (98) to the range 
0 < s < 1. 

Figure 4 gives a plausible curve for g(iq), and in addition the right 

Hco 

side of (96) as j drj . . . and g ' ( iq ) . The scale is selected in such 

Jo 

a manner that the two, last curves coincide for large tj. Considerable 
differences then exist for small tj but there the equation (96) also can 
no longer be correct. If one substitutes the function g(q) thus 
obtained into (91) and neglects the fluctuations , there results for k : 


K 



ds(-lg s)cp(s) 



&i}g(q)g(t]s) = °*98 


( 100 ) 


This crude estimation therefore gives the correct order of magnitude 
for k, but the exact value may well differ from O.98 by as much as a fac- 
tor 2. The calculations of this section thus have not led to an exact 
calculation of the constant. k but they did provide a qualitative mathe- 
matical representation of the processes on which the energy dissipation 
is based. Perhaps it will be possible to arrive at a rather exact experi- 
mental determination of k by means of a comprehensive discussion of the 
various experiments of Simmons, Dryden (cited before), Prandtl^-5, and 
others regarding the spectrum and the damping of turbulence . 

Translated by Mary L. Mahler - 

National Advisory Committee 
for Aeronautics 

•^prandtl, L. : Proc. VI. Intern. Congr. f. Appl. Mech. Cambridge, 

Mass., 1938; P- 340. 
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Figure 2. - The turbulent energy distribution as a function of the wave 

number. 







